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Abstract. The ionic Hubbard model on a cubic lattice is investigated using analytical approximations, the 
DMFT and Wilson's renormalization group for the charge excitation spectrum. Near the Mott insulating 
regime, where the Hubbard repulsion starts to dominate all energies, the formation of correlated bands is 
described. The corresponding partial spectral weights and local densities of states show the characteristic 
features, of a hybridized-band structure as appropriate for the regime at small U, which at half-filling is 
known as a band insulator. In particular, a narrow charge gap is obtained at half-filling, and the distribution 
of spectral quasi-particle weight reflects the fundamental hybridization mechanism of the model. 

PACS. 71.27.+a Strongly correlated electron systems; heavy fermions - 71.10.fd Lattice fermion models 
(Hubbard model, etc.) - 71.10.+h Metal-insulator transitions and other electronic transitions 



1 Introduction 

Investigations of the ionic Hubbard model have led to an 
ongoing debate about different types of phases with and 
without broken symmetries and transitions between them 
. Reliable information has so far been obtained mainly 
for the one-dimensional case, where Quantum Monte Carlo 
methods |H] or the density matrix renormalization group 
0] are applicable. From these it became clear, that a site 
dependent variation of atomic energies can induce critical 
behaviour and long-ranged correlations which are absent 
in the homogeneous Hubbard model. This simple version 
of the model in one dimension and for the particularly in- 
teresting case of half-filling is a correlated Mott insulator 
for all values U > of the local Coulomb repulsion, with 
a gap for charge excitations but none for spin excitations. 
Its critical bond-bond correlations apparently freeze out 
at low T in a parameter range U% < U < U2, when the 
model is chosen ionic. For U above U2 spin-spin corre- 
lations seem to dominate, so that in a regime, where U 
is sufficiently larger than the local energy variation, the 
same type of Mott insulating phase is approached as in 
the nonionic model. 

The homogeneous Hubbard model in higher dimen- 
sions, namely d = 3, likewise shows an interesting phase 
diagram, which e.g. was calculated in the framework of 
dynamical mean field theory [5]. Apart from regions possi- 
bly destabilized by phase separation, magnetic phases pre- 
vail for large values of U near half-filling, whereas upon 
doping the Mott insulator in form of a correlated para- 
magnet is increasingly stabilized. Studying the paramag- 
netic state can serve as a good starting point even in the 
regime, where the model develops magnetic order. It con- 



tains essential information about one particle excitations 
and their residual interactions. In particular, when the low 
lying quasi-particles form well defined bands, one may pos- 
sibly apply concepts known from the theory of weak band 
magnetism and from Stoner theory and thus describe the 
magnetic phase in terms of exchange splittings of heavy 
quasi-particle bands, at least near the phase transition line 
|fill4j . It is therefore desirable to investigate a correlated 
paramagnetic state for the Hubbard model at interme- 
diate U and small ioniticity, where formation of low ly- 
ing quasi-particle bands can be expected via an interplay 
of local correlations and inhomogeneous local energies. 
We will thus concentrate on one-particle excitations in a 
Fermi liquid phase and study the corresponding spectral 
functions via appropriate many-body techniques. Work on 
possible magnetic instabilities in the correlated bands is 
in progress. 

In the regime near the metal insulator transition, where 
U becomes comparable with the bandwidth as derived 
from nearest-neighbour transfer, the Hubbard model near 
half-filling shares an important feature with the Anderson 
lattice model: Bands of heavy quasiparticles with long life- 
times form at low temperatures. Although in both cases 
this process is driven by the strong local interactions and 
involves complicated many-body correlations, some fea- 
tures of the quasi-particle band structure seem to be linked 
to properties already inherent in the one-body terms of 
the Hamiltonian. This applies e.g. to the volume of the 
Fermi surface due to Luttingers theorem in the Fermi liq- 
uid state, but even more detailed structure may be pre- 
served as can be inferred from existing work on these mod- 
els |6I7I9I11| . Particularly intriguing is the case of the spec- 
tral composition of quasi-particles at different wavenum- 
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bers, which one could expect to be largely determined by 
the structure and symmetry of the fundamental hybridiza- 
tion terms of these models. We will show that indeed this 
also occurs in the ionic Hubbard model. 

The motivation to link a study of the ionic Hubbard 
model to results for the Anderson lattice is twofold: As 
will be discussed in the following, ionicity of the Hub- 
bard model causes similarities concerning the Brillouin 
zone and the band structure. Furthermore, new features 
to be expected due to these similarities require a certain 
standard of many-body techniques, a lecture which was 
originally learned in connection with the Anderson lattice. 
Application of the Friedel sum rule and of the Luttinger 
theorem to the Anderson lattice model sheds some 
light on the density of states (DOS) in the region of low 
lying quasi-particles |13| . In particular for one conduction 
band and two electrons per site the noninteracting version 
of the model has a hybridization gap at the Fermi energy, 
and the Fermi surface coincides with the boundary of the 
first Brillouin zone, a situation which should be preserved 
when interactions are included. More generally, one may 
argue that a coherent action of the Kondo effect should 
lead to a gap at the energetic position where the many 
body resonances form on each lattice site. Due to Friedels 
sum rule it will depend on the localized charge whether 
this position is near the Fermi energy and thus leads to 
a charge excitation gap. Indeed, many-body calculations 
using various techniques have demonstrated that a pic- 
ture of hybridized quasi- particle bands is an appropriate 
description of the low temperature regime of this model 
|14ll5llrill7ll8| . 

The Hubbard model likewise contains interactions and 
hybridization terms, the latter usually being addressed as 
transfers to nearest neighbours. The noninteracting ver- 
sion gives one simple tight binding conduction band, and 
the generic case puts the Fermi level somewhere near the 
band center. Correspondingly, the Fermi surface lies well 
inside the first Brillouin zone and a gap is not expected in 
the interacting case, too. The situation can change drasti- 
cally when an ionic Hubbard model is considered. For a bi- 
partite, simple cubic lattice with local one particle energies 
differing between the two sublattices, the first Brillouin 
zone is halfed, and with one electron per site the Fermi sur- 
face coincides with the new, smaller zone boundary. Then 
it can again be expected that the symmetry-breaking ef- 
fective field produces enough Bragg scattering to lift the 
degeneracy at the zone boundary and thus leads to a gap 
like in the Anderson lattice with two electrons per site. We 
would therefore expect on the basis of this similarity that 
at least near half-filling quasi-particle bands form at low 
T, which reveal the underlying hybridization mechanism. 
In addition, it would be interesting to find out, whether 
the wave vector dependence of the corresponding matrix 
elements is preserved in the structure of quasi-particles. 

This expectation, in fact, is supported by a theoreti- 
cal approach, which has proven extremely useful for such 
models with dominant local interactions. The approach 
rests upon the picture of effective sites |T3|. which react 
to a surrounding medium which in turn is formed by these 



same local objects. Links between them are established via 
the same elementary transfer processes which lead to band 
formation via tight binding of noninteracting electrons. 
Scattering of quasi-particles by the effective sites thus re- 
flects in lowest order just the k-dependence of these hy- 
bridization matrix elements, whereas the important part 
of the selfenergy remains local, i.e. k-independent. Quite 
a long line of approximations have been formulated along 
these lines, from the ATA ^S] an d the Renormalized Band 
Theory to the XNCA [T7j and the Dynamical Mean 
field Theory (DMFT) [T^J, which has become very popu- 
lar. It possesses a sound justification and has formed the 
basis of extensive studies of Hubbard and Anderson mod- 
els [0] . Generally, quasi-particle bands of the type outlined 
above are derived with these methods supporting our view 
and furnishing a good prospect for generalizations e.g. to 
the ionic Hubbard model. 

Early approaches to the Anderson lattice model like 
the ATA and the LNCA used simplified versions of the 
Non-Crossing Approximation (NCA) at finite U for the 
calculation of scattering processes by the effective sites. 
They clearly pointed to a quasi-particle band structure 
with a hybridization gap. A more consistent description of 
the effective medium, as contained in the XNCA, put this 
in doubt for a while, because subtle cancellations between 
local and nonlocal contributions to the selfenergy do not 
occur in the right way due to shortcomings of the NCA. 
When combined with a better local method the selfconsis- 
tency cycle of the XNCA produces the gap, too. In com- 
bination e.g. with the Numerical Renormaliztion Group 
(NRG) ^1 this has become a reliable tool, at least for 
low temperatures and excitation energies. The simplified 
NCA, as an analytical tool, can be expected to become 
still more useful here, when an extended or full version is 
manageable in numerical calculations. In connection with 
the Hubbard model this conceptual frame is known as the 
Dynamical Mean Field Theory (DMFT) and will be used 
for the essential calculations in section 3 of this paper. 
Since less ambitious approximations in the spirit of Hub- 
bard I or with analytical impurity solvers like SNCA in the 
DMFT-cycle are helpful to understand the initial stages of 
correlation effects and to span Luttinger's szenario from 
the noninteracting to the fully interacting case, we first 
give some elementary calculations demonstrating rough 
overall features of the expected hybridized band structure. 
The reader mainly interested in the correct final form of 
the results should proceed to the second half of section 3. 

In the following section 2 we will introduce our model 
Hamiltonian and discuss shortly the techniques used to 
solve the model. The transition from the energy levels of 
isolated ions towards the fully developed picture of corre- 
lated quasi-particle bands will be performed in two steps: 
At first, it is instructive to study the splitted bands fur- 
nished by a generalized Hubbard I-approximation, called 
Free Theory in the following for brevity |33| . which ne- 
glects all two particle correlations except for the purely 
local ones. It allows for a discussion of partial densities 
of states (DOS), interaction - and hybridization - induced 
band splittings and van Hove singularities in the absence 
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of lifetime effects. Section 3 proceeds to a generalized 
DMFT-scheme, which includes the complete set of local 
correlations and is self-consistent on the lattice. Calcula- 
tions are carried out with an analytical impurity solver, 
a simplified version of the finite U - NCA '_'() , and al- 
ternatively with the NRG. Whereas the former allows for 
a qualitative study of lifetime effects and is particularly 
useful at higher excitation energies, a clear cut picture of 
hybridization effects in the region of the many-body reso- 
nance is only achieved with the latter. The final section 4 
contains a short discussion and concluding remarks about 
future perspectives. 



2 The bipartite Hubbard model and 
correlated bands 

We build our model with two sorts of s-shells with one 
particle energies £a,£b — £a + A and local Colomb ma- 
trix elements Ua,Ub- Each of them is placed onto one 
sublattice La, Lb of a three-dimensional simple cubic lat- 
tice, so that the nearest neighbours of one sort belong to 
the other. The Hamiltonian is: 
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We have introduced the local occupation operators 
nAja = a ^ a a j<? an d n B£a = b^ a ^>ia for sites on each sublat- 
tice and absorbed the (common) chemical potential into 
the one particle energies. A scheme of the resulting ener- 
gies of local shells is shown in Fig. which also gives an 
impression of the energy regime we have in mind when the 
indicated position of the Fermi level and the zeroth order 
bandwidths are recognized. 

This atomic limit is described by the Hamiltonian Ha+ 
Hb ■ If the transfer term is added instead of the interac- 
tions, i.e. when 



H (o) = H (0) 



H 



(0) 



Hab 



(2) 



is considered, a picture of hybridized tight binding bands 
for noninteracting electrons emerges, which is easily vi- 
sualized: Assuming first A = 0, the single tight binding 
band of a homogeneous (e A = £b) lattice is folded back 
corresponding to a halfing of the Brillouin zone. The ef- 
fect of A then shifts the remaining original and backfolded 
pieces apart, producting gaps due to Bragg-scattering be- 
tween the boundaries of the new reduced zone. Fig. Ogives 
an impression of this band structure, using a single co- 
sine for the dispersion function er' — 2icos(fca) (a = 
lattice constant, t < 0) as appropriate for one dimension. 



Fig. 1. Scheme of the local shells and the energy regime used 
for the bipartite Hubbard model. The expected bandwidth is 
indicated by the vertical delimeters. The Fermi level lies well 
separated between the lower and upper bands. 

Shown are the spikes corresponding to the partial DOS 
a + (ui + i5) (pb analogous) 

' kcr 

(3) 

derived separately for excitations on each sublattice. 

It is easy to see (e.g. perturbatively with respect to t 
or with the equation of motion method) that the Greens 
functions appearing here are given by 



p A (k,u) = Im G 

IT 



kfr 

b + analogous) 



(4) 



where the counterparts for the noninteracting atomic limit 
have to be inserted, i.e. 



Gf{z) = [z-e A ] 



{G^g analogous). (5) 



Whereas the added spectral weight p = pa + Pb would 
simply exhibit two (5-spikes with weight h at the energies 
of the upper and lower band for each fixed value of crystal 
momentum k and spin a, the distribution of the partial 
spectral weights is more interesting, following 

PA ^w) = (^(k) • 8(w - 4 +) ) + Ci +) (k) • 8(w - e^), 

A 
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_(±) 



(.A - 2 ± 



2V^ 2 + (24 0) ) 2 



£A 



eB±\/A2 + (2e i ° ) y 



(6) 



4 



Torben Jabben et al.: Charge gaps and quasiparticle bands of the ionic Hubbard model 




2a €A + 2i 



Fig. 2. Tight-binding approximation of the one dimensional 
partial DOS for the A and B sublattice. 

For pb one has to exchange Cg = Ca> Cb = Ca + ^- ^ e 
gap at the zone boundary equals that of the local states, 

i.g. e<+) - e(-) = \/z\ 2 + (24 0) ) 2 =A for k = ±£, and 
at these two points the bands terminate in the local en- 
ergies — £ B and 4 = e A for fc = if-, inducing 
complete repulsion of admixture from the corresponding 
other state. Therefore the upper band has full S-weight at 
the zone boundaries and less in the middle, whereas the A 
weight goes to zero at the boundaries and grows towards 
the middle of the Brillouin zone and the other way round 
for the lower band. This simply reflects the shape of the 
dispersion function which is zero at the boundaries 
and has maximal absolute value \2t\ in the zone center, 
since this function also gives the effective hybridization of 
electrons with momentum k in transfer processes between 
a site and all of its nearest neighbours. 

It is worthwhile to collect these simple facts, since the 
effects are relevant also for the more refined band structure 
derived with better approximations. A generalization of 
the Hubbard I-approximation, named Free Theory, can 
be obtained by simply substituting the Greens functions 
of the isolated interacting local shells, i.e. 

r (2) 

G<-}>(z) = -^— + ^ TT (g£ } analogous), 

z - sa z - sa - Ua 

(7) 



for the G ( X> B given by Eq. (5) in the explicit expressions 
(4). The resulting bands are derived from a polynomial 
of fourth order, the zeroes of which are easily determined 
numerically, as well as the corresponding partial spectral 
weights, too. Here we give the explicit dispersion for the 
special case Ua = Ub = U and T = where Q 1 ' = 
Ca ' 1 = Cb = Cb — \ m tne situation depicted in Fig. ^ 

e M = l -[e A + £B + U± 

\Ja 2 + U 2 + (2e^ 0) ) 2 ± 2 y/ {AUf + (Us^) 2 + (e^ 0) ) 4 ) , 
with 

m = 1=(+, +), m = 2=(+, -), 

m = 3=(-,-),m = 4=(- +). (8) 

The band structure is shown in Fig. 01 where the spectral 
weights are again resolved in contribution from sublattices 
A and B. Throughout the paper, all energies will be given 
in units of the hopping parameter \t\. 



PB(k,U>) 



0.6 r 




Fig. 3. Partial one dimensional DOS of the A and the B sub- 
lattice from Free Theory calculations for €a = —2, es = —1.66, 
Ua — Ub =4. Due to the energy difference A — es — sa each 
Hubbard band splits in two subbands with hybridization gaps 
of width A. 
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One recognizes the same distribution of weight as in 
the noninteracting case before. The important new fea- 
ture occuring in the Free Theory is the interaction-induced 
additional splitting of the hybridized bandstructure, for- 
merly seen in Fig. [21 which is in complete agreement with 
the phenomenon known from Hubbards own study: Each 
of the bands splits in two subbands of approximately half 
the spectral weight, and differing in energy by roughly the 
Coulomb repulsion. Interesting is also the corresponding 
local DOS, which is obtained via summation over all k. A 
glance at the upcoming Fig.0]shows typical van Hove sin- 
gularities at the bandedges. In the present bipartite case 
they may always be attributed to one of the sublattices 
only, since for the other sublattices the spectral weight 
vanishes at the corresponding energies. Concluding this 
section, we comment shortly on the validity of Eq. (4) and 
similar generalizations of the Hubbard I-approximation. 
The essential point is the neglect of all nonlocal correla- 
tions and thus the restriction to the most important lo- 
cal interaction effect. This is contained in the fractional 
form of the Greens function (7), which is easily general- 
ized to a more complicated local level structure. Using 
this Greens function, the approximation simply proceeds 
as in the noninteracting case, involves the Wick theorem 
and all benefits thereof. This may also be obtained from 
a decoupling in equations of motion. 



3 Selfconistent approximations for effective 
sites and low lying quasi-particle states 

Corrections to the Free Theory can be organized perturba- 
tively with help of local cumulant interactions serving as 
vertices. In order to tackle the infrared problem connected 
with degeneracies in the metallic regime resummations to 
infinite order are necessary |21I22| . The well established 
strategy of DMFT uses the concept of effective sites, self- 
consistcntly coupled to an environment, which is built up 
by transfering electrons between these same objects. This 
is accomplished by using an appropriately highly devel- 
opped technique for an impurity problem and by feeding 
it with effective propagators for electrons leaving and en- 
tering and being scattered in between by other effective 
sites of the lattice. It is well known how this concept is for- 
mally implemented in an elegant manner as a small set of 
equations |17I18| : we will just state here its generalisation 
to the present situation with a few comments concerning 
the role of the quantities appearing. 

The two types of lattice sites experience different forms 
of local one-particle excitation processes, which as before 
are described via a Greens function in the form 

G A {z) = G a ,^ a + a (z) = \ ^G aka a +^{z) (G B analogous). 

k 

(9) 

For a description of the propagation process inherent in 
the k-dependent Greens function under the sum a divi- 
sion into free transfer processes and effective local parts 
is helpful. This enables one to apply the same reasoning 



as was used to derive Eq. (4) and thus gives an analogous 
result here: 

G a k „.a+ M = [Ga(z)- 1 (e k ) 2 G s (z)l _1 

(G bk a + analogous). (10) 

The effective local Greens functions Ga and Gb contain 
corrections for processes, where particles leave and en- 
ter the site along an irreducible loop of transfers. Since 
such loops on the one hand contribute to the pseudolocal 
Greens functions Ga and Gb and on the other are pro- 
duced by the free random walk through the lattice lead- 
ing to the result Eq. (10), these corrections are necessary 
to avoid overcounting. They are collected in quantities 
Ta(z) and Tb{z), which enter in the following way: 

Ga(z)^ 1 — Ga(z)^ 1 + Ta{z) (Gb(z) analogous), 

(11) 

where Ga results from Eq. (9). In this way a closed self- 
consistency cycle is obtained, in the present form with 
coupled A- and £?-sublattices, which has to be supple- 
mented as mentioned by an impurity theory for Ga and 
G B - In the framework of the DMFT |18l8l9j . the quanti- 
ties Ta/b( z ) are usually called the dynamical Weiss fields 
of the theory. Recently it was pointed out fTI]] that they 
determine how much the effective local inverse Greens 
function Ga/b^z) -1 has to be deformed to reproduce the 
k-summed lattice Greens function Ga/b ■ 

Analytical impurity solvers have been derived from the 
noncrossing approximation (NCA) of the Anderson impu- 
rity model |2.'ij . Their virtues and shortcomings are well 
known [21] . A full version of the NCA for finite Coulomb 
repulsion has been presented long ago |2()j , but for a quick 
orientation or for complicated applications usually a sim- 
plified version, the SNCA is employed, which needs much 
less numerical effort. Also improvements of the full NCA 
have been presented |25I26| , which however consume even 
more time and effort. The SNCA, and to a lesser degree 
the full NCA, loose their relliability in the regime of very 
low temperatures and excitation energies. They describe, 
however, rather well the overall structure and can be used 
down to temperatures and excitation energies of the or- 
der of the many-body scale, connected with the infrared 
problem, and somewhat below. In particular the full NCA 
reproduces this nonperturbative scale with quantitative 
accuracy 120] • We have employed here the SNCA in order 
to obtain a first test of the selfconsistency cycle outlined 
in Eqs. (9) to (11). Thereby, Ga and Gb are calculated us- 
ing the loops Ta and Tb as input, afterwards one obtains 
Ga and Gb using (11) with these Ga and Ta, likewise 
for Gb , and finally a new Ga and Gb is produced with 
(10) inserted into (9). These Ga and Gb give rise to new 
loop-expressions Ta and Tb, again via Eq. (11), which in 
turn allow for a new impurity calculation of Ga and Gb- 

In Fig. ^ results of this iterative procedure are shown 
and compared to the corresponding quantity calculated 
in the Free Theory Fig. fjjc) contains the complete local 
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DOS; 

Pab{u) = -—Im[G A (v + iS) + G B (w + iS)] = 

= ±\p A (w)+PB(w)]. (12) 

and compares it to the corresponding quantity calculated 
in the Free Theory. 




-4 -2 2 4 

OJ 

(C) 

Fig. 4. Comparison of the three dimensional local DOS 
obtained from Free Theory (solid line) and DMFT+SNCA 
(dashed line) for <l a = -2.07, e B = -1.89, Ua = U B = 3.96, 
= 20. 

It is surprising, how much lifetime effects, i.e. the broad- 
ening of DOS peaks at fixed k, smear out the prominent 
features of the DOS, which have been discussed in Section 
2. The hybridization gaps now only show up as shallow 
dips and the van Hove singularities have disappeared. The 
prospect, however, to obtain hybridization-induced gaps 
in the narrow quasi-particle peak around the Fermi level 
oj = are unbroken: The Fermi liquid state guarantees di- 
vergent lifetimes for T, u) — > 0. It seems thus encouraging, 
that the calculation indeed shows some precursor of split- 
ting near lu = 0, which is more pronounced in the partial 
DOS shown in Figs. Eta) and Efb). Unfortunately, our 
selfconcistency cycle with the simplified SNCA tends to 
become unstable in the most interesting region and does 



not allow for a more precise investigation of this interest- 
ing effect. One should notice from an inspection of Fig. EI 
however, that some of the characteristic features of the bi- 
partite lattice, which have been discussed before, survive 
even in such a locally complete calculation, although only 
in a smoothed form. 

In order to answer reliably the question for hybridiza- 
tion structure in the low energy quasi-particle domain we 
have finally employed Wilsons renormalization group as 
the impurity solver in the selfconsistency cycle |27I28I3()] 
defined above. In the meantime, this has become a stan- 
dard procedure, which essentially improves on the low 
energy-low temperature results, but due to its logarithmic 
partition of excitation energies around uj = tends to sup- 
press important DOS structure elsewhere This latter 
statement is substantiated in Fig. 5, viewing the smooth 
and featureless main resonances near ea and ea + Ua, 
likewise for B. At the Fermi level ui = 0, on the other 
hand, this T = calculation reveals a clear and complete 
hybridization gap in the quasi-particle DOS at the low- 
est excitation energies. Although the present calculation 
uses parameter values Ua = Ub = 3.96, ea — — 2.07,es = 
— 1.89 near the Mott transition, we expect at least for half- 
filling this gap to exist at the Fermi level for all values of 
U in the Fermi-liquid state, as we have motivated above 
with reference to the Luttinger theorem. 
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Fig. 5. Three dimensional DOS from DMFT+NRG for e A = 
-2.07, e s = -1.89, U A = U B = 3.96, T — > 0. The NRG 
treatment of the ionic Hubbard model shows a clear quasi- 
particle hybridization gap at the Fermi level. For the NRG 
we used A — 2.3 as discretization parameter |27| and kept 
N B t = 700 states for each NRG iteration. 



Torben Jabben et al.: Charge gaps and quasiparticle bands of the ionic Hubbard model 



7 



PAB(k, W) 
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Fig. 6. [Ill] k resolved DMFT+NRG spectra with the same 
parameters as in Fig. 03 The top picture shows the whole spec- 
trum for the combined DOS while the other two picures show 
the partial DOS for the A and B sublattice for the region 
around the quasiparticle bands at the fermi level. 



We have suggested above, that the quasi-particle states 
are shaped by hybridization in a similar way as the 
Hubbard-split bands of width W oc t in the high energy 
region. It would therefore be consistent, if the size of this 
gap compared to the one seen e.g. in Fig. 3 is scaled 
down roughly by the same factor of 7 = T* jW as the 
whole band structure. T* is a characteristic low energy 
scale caused by the strong correlations, which should be 
roughly of the same magnitude as the enhanced Kondo 
scale of the nonionic version of the model. As was pointed 
out in section 2 the gap in Free Theory equals A, so that 
a quasi-particle gap size A* of order 7 ■ A can be ex- 
pected. Whereas this order of magnitude actually is found 
in our calculation, and A* clearly vanishes when A goes 
to zero and increases with growing A, a strict proportion- 
ality is not to be expected. Lifetime effects in the DMFT 
calculation start to smear features in the DOS when the 
distance to the Fermi level becomes larger and tend to 
narrow a gap situated there. In addition, the existence of 
a quasi-particle gap causes a change in the scale T* as 



compared with the nonionic case. It would, however, be 
interesting to study the variation of A* near the metal- 
insulator transition or away from half-filling, although one 
has to expect numerical difficulties. A further similarity 
concerns the distribution of spectral weight in the partial 
DOS, which formerly was discussed in connection with the 
structures at higher energies. Even without k-resolution 
the spectral weight for the A-lattice near 10 = in Figs. 
5(a) and 5(b) resembles the asymmetric shape observed 
near the hybridization gaps at higher excitation energy 
in Figs. 4(a) and 4(b), and also the k-dependence reveals 
strong similarities. This is demonstrated in Fig. 6, where 
the upper part gives an overall view of the total weight 
Pab (k, u>) and the two lower parts magnify the region very 
near to the Fermi level oj = 0. 

Regarding the regions at large excitation energies ui 
in the upper figure, considerable broadening and spread- 
ing of the resonances is observed compared to those of 
the Free Theory, see e.g. Fig. 3. Moreover, band splitting 
due to the ionic field A = eb — s A is not visible any- 
more and it seems that parts of the broadened band struc- 
ture even have disappeared. Although the NRG treatment 
looses much accuracy away from the Fermi level due to the 
logarithmic discretization of energies, we attribute these 
findings to the dominance of scattering processes in this 
region due to the blocking effect on the effective sites. 
For the ionic version of the Hubbard model this seems 
to be of particular importance. Without the ionic field A 
one may define the positions Wk of broadened bands as 
usual via the disappearance of the real part in the de- 
nominator of one particle Greens functions, i.e. solving for 
3?e G(cj k + ^)" 1 -e k = w k - (e k + 2 + sft e £(uj k )j = in 

the Hubbard model, where we have separated a Hartree- 
part y from the selfenergy of the effective Greens function 
G(z) and explicitly used its k-independence in the effec- 
tive site picture. Although 3?e E(oj + id) bears a strong 
w-dependence, generally three solutions appear for the 
Hubbard model in the Fermi liquid regime which rep- 
resent the two original [/-split bands and the low- lying 
quasi-particle band. Level broadening occurs separately 
via Im£(uik + iS). For the ionic Hubbard model the de- 
nominator of the Greens functions JU| mixes contribu- 
tions of the two kinds of effective sites in a way that 
stresses the importance of scattering even more, i.e. one 
has to solve 

(u} k - e A - S A (uJw)j \ Uk - £b - ^ - 5Rei7_B(w k )^ 
— /m^i(wk + iS) ■ ImSsi^k + ifi) — e k = 0. (13) 

Roughly speaking, ImEA{u-k+i$) ~ JmlTB^k+i^) rs As 
is a large energy shift for w k of order ea,£a + U , so that 
only the solutions, which are maximally removed from the 
Fermi level tend to remain in the high energy region after 
inclusion of the imaginary part of the self energies. In the 
low energy region, on the other hand, this effect is small- 
est near the Fermi level, which around half-filling makes 
the zone boundaries favourable for the existence of well 
defined quasi-particle bands. Thus we observe more pro- 
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nounced resonances there as shown in the two lower parts 
of Fig. 6. 

Apparently, the dominant parts of the quasi-particle 
bands reproduce the distribution of main spectral weight 
on the respective sublattice. The parts on the opposite 
side of the Fermi level, on the other hand, are shaped by 
admixtures, e.g. the spectral weight at the zone bound- 
aries does not vanish anymore. Near the position of the 
ionic energies the NRG treatment combined with the ef- 
fects described before leads to a less detailed band struc- 
ture. This region is presumably described in a somewhat 
better way by the NCA calculation, see e.g. Fig. 7. Thus, 
our expectations that fundamental hybridization mecha- 
nisms, contained in the model, bear a relevance on the low 
lying states with many-body character, is supported. We 
have compiled SNCA and NRG results in our last Fig. 7 
in order to provide an overall picture of the local DOS 
and both of its sublattice contributions. This also serves 
to underline our remarks about the usefulness of both cal- 
culational schemes in conjunction. 




4 Conclusion 

We have calculated correlated band structures of the ionic 
Hubbard model for a region of intermediate local repul- 
sion at and near half-filling. The quasi-particle bands of 
low lying one-particle excitations show clear signatures of 
the underlying transfer (" hybridization" ) mechanism con- 
tained in the noninteracting part of the Hamiltonian, in 
particular with respect to a charge excitation gap and the 
distribution of spectral weight. We have motivated and 
interpreted these findings with reference to the Luttinger 
sum rule and to the dominantly local nature of the selfcn- 
ergy and have drawn some parallels to gap formation in 
the Anderson lattice. Applying successively better approx- 
imations we tried to elucidate the transformation from a 
local ionic picture to essentially itinerant quasi-particle 
degrees of freedom and the correspondance to an interme- 
diate tight binding scenario. 

Using the SNCA in the DMFT sclfconsisting cycle gave 
a good overall impression of the band structure and al- 
lowed for some qualitative conclusions about the low lying 
excitations. A satisfactory picture of the low energy region 
at low temperatures was nevertheless only achieved by ap- 
plication of the numerical renormalization group. Whether 
a full or improved NCA would allow for a more analytical 
approach and would lead to sensible results in this region, 
too, is still to be demonstrated. Altogether, the scope of 
controlling and understanding such calculations for com- 
plicated systems with more ionic species and realistic p- or 
d-shells, in the local approach or e.g. in the LDA-DMFT 
scheme jllj . has been widened and improved: Some fea- 
tures of heavy quasi-particlcs seem definitely closely con- 
nected to a noninteracting system. Since the whole field 
at present rapidly develops to a stage, which is technically 
very much involved, this insight should be useful for fu- 
ture calculations. 

As far as the existence of additional long ranged corre- 
lations and corresponding order parameters are concerned, 
a possible instability of the Fermi-liquid state considered 
here towards magnetic order or towards phase separation 
should be investigated. Like in the homogeneous Hubbard 
model [2] knowledge of the one particle excitations is 
an important prerequisite for this program. Work on two 
particle properties within the DMFT approach, and in 
particular on magnetic and charge susceptibilities [■V2\ . 
is in progress. Without further calculations our results 
also demonstrate an asymmetry of occupation numbers 
n>A - "b / away from half-filling as a consequence 
of the asymmetric distribution of spectral densities, see 
e.g. Fig. 4(a, b). 



Fig. 7. Comparisons DOS obtained from DMFT with NRG 
and NCA for a 3d simple cubic lattice. The insets show the 
region around the Fermi level in greater detail. The parameters 
used are the same as in Fig.QO for the NCA a finite temperature 
1/T = /? = 20 was used. The NCA resolves much better the 
high energy features of the spectra but is not able to clearly 
resolve the quasiparticle hybridization gap as the NRG. 
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